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Abstract 

Hu, Paz and Zhang [ B.L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. D 45 

(1992) 2843] have derived an exact master equation for quantum Brownian 

motion in a general environment via path integral techniques. Their master 

equation provides a very useful tool to study the decoherence of a quantum 

system due to the interaction with its environment. In this paper, we give an 

alternative and elementary derivation of the Hu-Paz-Zhang master equation, 

which involves tracing the evolution equation for the Wigner function. We 

also discuss the master equation in some special cases. 
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I. INTRODUCTION 



Quantum Brownian motion (QBM) models provide a paradigm of open quantum systems 
that has been very useful in quantum measurement theory [1], quantum optics [2] and 
decoherence [3-5]. One of the advantages of the QBM models is that they are reasonably 
simple, yet sufficiently complex to manifest many important features of realistic physical 
processes. 

Central to the study of QBM is the master equation for the reduced density operator of 
the Brownian particle, derived by tracing out the environment in the evolution equation for 
the combined system plus environment. A variety of such derivation have been given [6-9]. 
The most general is that of Hu, Paz and Zhang [10,11], who used path integral techniques 
and in particular, the Feynman- Vernon influence functional. 

The purpose of this paper is to provide an alternative and elementary derivation of the 
Hu-Paz-Zhang master equation for QBM, by tracing the evolution equation for the Wigner 
function of the whole system. 



II. MASTER EQUATION FOR QUANTUM BROWNIAN MOTION 

The system we considered is a harmonic oscillator with mass M and bare frequency Q, 
in interaction with a thermal bath consisting of a set of harmonic oscillators with mass m n 
and natural frequency uj n . The Hamiltonian of the system plus environment is given by 

" = w + \ Ma ^ 2 + ? (4 + \ m ^A + "? c " 5 " • (1) 

where q,p and q n ,p n are the coordinates and momenta of the Brownian particle and oscil- 
lators, respectively, and C n are coupling constants. 

The state of the combined system (1) is most completely described by a density matrix 
p(q, Qi'i 4 1 q'ii t) where g» denotes (qi, ...(/at), and p evolves according to 

P=~[H,p\. (2) 
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The state of the Brownian particle is described the reduced density matrix, defined by 
tracing over the environment, 

p r (q, q', t)= Y[(dq n dq' n 5(q n - q' n ))p{q, qu q', q-, t) . (3) 

J n 

The equation of time evolution for the reduced density matrix is called the master equation. 
For a general environment, Hu, Paz, and Zhang [10] derived the following master equation 
by using path integral techniques: 

dp r h 2 (d 2 p r d 2 p r \ l,, n2 / 2 a\ 

+^M5Q 2 (t)(q 2 -q' 2 )p r 
-iMV{t)h{t){q-q'fp r 

+ ar (t )/( t )(,-,o(^ + ^). (4) 

The explicit form of the coefficients of the above equation will be given later on. The coef- 
ficient 5Q 2 (t) is the frequency shift term, the coefficients T(t) is the "quantum dissipative" 
term, and the coefficients T(t)h(t), T(t)f(t) are "quantum diffusion" terms. Generally, these 
coefficients are time dependent 

and of quite complicated behaviour. 

We find it convenient to use the Wigner function of the reduced density matrix, 

W{q,p,t) = ±fdueWfr(q-±q+±i) . (5) 
Taking the Wigner transform of (4), we obtainF] 

d(vW) d 2 W 
+ 2T(t)^-l + HMT(t)h(t) — 

d 2 W 

+ sr( ( )/ (i )_. (6) 



1 We believe Eq. (2.48) in Ref. [10] contains some incorrect numerical factors. 
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The inverse transformation of (5) is given by 

Pr (q,q',t) = J dpe-^-^W • (7) 

Our strategy for deriving the master equation (4) is to derive the Fokker-Planck type 
equation (6) from the Wigner equation for the total system. The master equation can be 
obtained from the Wigner equation for the system by using the transformation (7). 

We shall make the following two assumptions: 

(1) The system and the environment are initially uncorrelated, ie. the initial Wigner 
function factors 

Wo(q,p;qi,Pi) = W^q,p)W^q i ,p i ) , (8) 

where Wq and Wq are the Wigner functions of the system and the bath, respectively, at 
t = 0. 

(2) The heat bath is initially in a thermal equilibrium state at temperature T = (/cb/3) -1 . 
This means that the initial Wigner function of bath is of Gaussian form, 



w»=nw, 



b 

nO 



= J] N n exp{ — tanh(-Hu n f3)H n } , (9) 



where H n is the Hamiltonian of the n-th oscillator in the bath, 



Hn = + ^ . (10) 



In addition, one easily see that the initial moments of the bath are 



<<Zn(0)> = (p„(0)) = , (11) 

(<7n(0)<? m (0)> = (if m^n), (12) 

(Pn(O)Prn(O)) = (if 171 ^ Tl) , (13) 

(?n(0)p m (0) + p m (0)g n (0)> = , (14) 



and 



(?n(0)) = ^— coth(^o; n /3), 

Zlll n UJ n z 

(Pn(°)) = ^hm n u n coth(^fuv n p) . (15) 

For the QBM problem described by (1) and (2), the Wigner function of the combined 
system plus environment satisfies 

dW p dW „ - dW 

-TT- = - — + MQ 2 q—— 
dt M dq y dp 

^\m n dq n n n n dp n ) 
^„ f dW dW\ 

+?c „^_ +9 _j. (16) 

By integrating over the bath variables on the both sides of the above equation , one obtains 
dW p dW %/rra dW f -■— |- J <9W^ 

where W^(?,p) is the reduced Wigner function and it follows from Eq. (3) that 

W(q,p)= [[dqifipiWfanq^Pi) . (18) 

J-oo i 

This definition is equivalent to Eqs. (3) and (5). The first two terms on the right-hand side 
of the Eq. (17) give rise to the standard evolution equation of the system. The last term 
contains all the information about the behaviour of the system in the presence of interaction 
with environment. 

In what follows, we shall demonstrate that the quantity 

G(q,p) = Y^Cn I X\dq t dp iqn W (19) 

n J i 

appearing (differentiated with respect to p) in (17) can be expressed in terms of IT^and its 
derivatives. To this end, we first perform Fourier transform of G(q,p) 

G{k, k') = J dqdpe ikq+ik ' p G{q,p) 

= £ C n / dqdp 1] dqidpi q n e ik ^ W(q, p;q t ,p t ). (20) 



It is well known that q(t),p(t) and q n (t),p n (t) are related to the classical evolution of 
their initial values q(0),p(0) and <?n(0),p n (0) through a canonical transformation: 



z(t) = U(t)z(p) , 



(21) 



where 



<t) = (q(t),q 1 (t)...q N (ty,p(t), Pl (t),...p N (t)) . 

Since the Hamiltonian (1) is quadratic, the Eq. (16) has the same form as the classical 
Liouville equation, so the solution of Eq. (16) is of the form, 



W t {z) = W {U-\t)z) . 



(22) 



Changing the integration variables into their initial values by this canonical transformation, 
we obtain 



G(k, k') = f dq(0)dp(0)l[dqi{0)dpi(0) 
J i 

fq(0) + g P (0) + £(/„<Z„(0) + 9nPn(0)) 



X 



xcxp 



xexp 



ik [aq(0)+(3p(0)) + J2( (0)+b nPn (0)) 

\ n 

ik'M (dq{0) + Pp(0) + ]T(d n g n (0) + b nPn {0)) 



x^(g(0),p(0))<(? i (0) I p i (0) 



(23) 



Here the coefficients /, g, f n , g n , a, (3, a n , b n are time dependent. Their explicit values are not 
required. 

Similarly, the Fourier transform of the reduced Wigner function is 

W(k, k') = J dqdpe ikq+ik ' p W(q,p) 

dg(0)dp(0)n<ift(0)dPi(0) 

i 

(0)+b nPn (0)) 
ik'M (aq{0)+$p{0) + J2( {0)+b nPn (0)) 



xcxp 



xexp 

xW°(q(0),p(0))W^q t (0),p t (0)) . 



(24) 
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Now compare G(k, k') and W(k, k'). They differ by the terms linear in q(0), p(0), g n (0), Pn(0) 
in the preexponential factor in G(k,k'). Consider the factors f n q n (0) and <7 n Pn(0) in 
G(k,k'). Since they multiply Wq (?i(0), j?i(0)), and since W / o(5'i(0),pj(0)) is Gaussian in 
?n(0),p n (0), the terms /„g„(0)W / 6 and <? n p n (0)Wo may be replaced by terms of the form 
dWQ/dq n (0), dWQ/dp n (0) up to time dependent factors. An integration by parts then may 
be performed, and these factors are then effectively replaced by multiplicative factors of 

Similarly, the factors fq(0), gp(0) in the prefactor in G(k,k') may be replaced by 
d/dk,d/dk' (plus some more factors of k and k'). Hence, it is readily seen that G(k,k') 
is a linear combination of terms of the form k, k' , d/dk, d/dk' operating on W(k,k'), with 
time dependent coefficients. 

Inverting the Fourier transform, it follows that 

dW dW 

G = A(t)qW + B(t)pW + C(t)— + D(t)— . (25) 

for some coefficients A(t), B(t),C(t), D(t) to be determined. This result immediately leads 
to the general form Wigner equation : 

dW p dW „ „ dW . . dW 



dt M dq dp dp 



III. DETERMINATION OF THE COEFFICIENTS(GENERAL CASE) 

Having found the functional form of the Wigner equation (26) of the Brownian particle, 
the next step is to determine the coefficients in the equation. Undoubtedly, there is more 
than one way to do this. Here we shall choose a way which is both mathematically simple 
and physically heuristic. Towards this direction, let us consider the time evolution of the 
expectation values of the system variables: q,p, q 2 ,p 2 and \{pq + qp). 

By using Eq. (16), we have 
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-^> = -Mfi 2 (g>-]TC< n <^>, (28) 

f t (q 2 ) = ^(pq + qp), (29) 

^(p 2 ) = -Mfi 2 (pg + qp) -2j2C n (pq n ) , (30) 

j t (pq + qp) = ^(P 2 ) - 2Mfi 2 (g 2 ) -2^C n (qq n ) . (31) 
Similarly, using Eq. (26) yields 

> = (32) 

±(p) = -(Mn 2 + A)(q)-B(p), (33) 

j t (q 2 ) = jj(pq + qp), (34) 

-(p 2 ) = -(Mfi 2 + A)(pq + qp) - 2B(p 2 ) + 2D , (35) 

j f (pq + qp) = ^(P 2 ) - 2(Mfi 2 + A)(q 2 ) - B(pq + qp) + 2C . (36) 

Since the evolution equations of the expectation values are confined to the system variables, 
the above two sets of equations must be identical. 
Now by comparing (28) with (33) we see that 

J2C n (q n )=A(q) + B(p). (37) 

n 

Similarly, by comparing (31) with (36), (30) with (35), respectively, we get 

£ C n (qq n ) = A(q 2 ) + ^( qp + pq )-C, (38) 

E C n(pqn) = ^{pq + qp) + B(p 2 ) - D . (39) 

The coefficients A,B,C,D may now be determined from (37)-(39) by regarding the 

expectation values (q), (q n q) etc. as expectation values of Heisenberg picture operators, and 
by solving the operator equation of motion. For simplicity, we still use ordinary notation to 
represent an operator without adding a hat on it. 
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The solution to the equation of motion may be written, 

q n (t) = a n q(t) + j3 n p{t) + ^2(a nm q m (0) + b nm p m (0)) , (40) 

m 

q(t) = aq (0) + (3p(0) + J> n g n (0) + b nPn (0)) , (41) 

n 

for some time-dependent coefficients a n , f3 n ,a nm ,b nm ,a, f3,a n ,b n . Note that q n (t) has been 
expressed in terms of the final, not initial values of q,p. By substituting Eq. (40) into Eq. 
(37), keeping (11) in mind, and comparing the two sides of the resulting equation, we have 

A = Y,C n a n , B = Y,C n [3 n . (42) 

n n 

Similarly, substituting (40) and (41) into (38) and (39), respectively, we get 

C = - E C n (a nm a m (q 2 m (0)) + b nm b m {p 2 m (0))) , (43) 

ran 

D = -Mj2C n (a nm a m (q 2 m (0)} + b nm b m (p 2 m (0)}) . (44) 

ran 

Here we have made use of p = Mq. The coefficients A, B, C, D are therefore completely 
determined by solving the equation of motion. We now do this explicitly. 
We have 

qXt)+{l 2 q(t) = -±-Y, C nqn(t), (45) 

n 

q n (t)+u 2 n q n (t) = -^q(t). (46) 
ui n 

The solution to Eq. (46) is as follows: 

/ , \ / n x , ,n . Pn(0)sin(w n t) 

9n(t) = 9n(0) COS{U n t) + 

-C n r d s SinMt ' s)] q -^. (47) 
Jo u n m n 

Combining (45) and (47) gives 

qit) + Q 2 q(t) + A jf dn^f - S )g( S ) = M , ( 48 ) 



where 

p n (0) sm(u n t) 



f(t) = ~J2 C n U„(0)cos(w n t) + 

n \ 

The kernel 77(5) is defined as 



(49) 



V(s) = ^(s) , (50) 

where 

r+°° 

7(5)= / dco cos(cjs) . (51) 

Jo CO 

Here J(w) is the spectral density of the environment: 

7H=£%-^)^L. (52) 

In order to get the expressions (40) and (41), we solve equation (48) with the following 
two different initial conditions: 

q(s = 0) = 9(0) ,q(s = 0) = ^ . (53) 

and 

?(* = *) = 9(0 , <K« = t) = ^ • (54) 

where t is any given time point. In doing so, we consider the elementary functions Ui(s)(i = 
1, 2) introduced by Hu, Paz and Zhang [10] which satisfy the following homogeneous integro- 
differential equation 

E(s) + fi 2 S(s) + — J° d\r]{s - A)S(A) = (55) 
with the boundary conditions: 

Ul (s = 0) = 1 , ui(s = t) = , (56) 

and 
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u 2 {s = 0) = , u 2 (s = t) = l. (57) 

The solution to equation (55) with the initial condition (53) is obtained as the linear 
combination of u±,u 2 , 

•w-(-.w-S-w)«w)^^- m 

The solution to equation (48) with the homogeneous initial conditions can be formally 
written as 



wis) 



± J o S drG 1 (s,r)f(r). (59) 
Where Gi(si, s 2 ) is the Green function which can be constructed in terms of Ui(i = 1, 2): 



s 2 ) = < 



m(32)M2(si)-M2(s2)m(ai) > 
«1(S2)"2(S2)-Wl(s2)"2(s2) 1 2 



otherwise. 
Note that s 2 ) as the function of si satisfies equation (55) with 

Gi(si = s 2 ,s 2 ) = 0, -^—Gi(si = s 2 , s 2 ) = 1 . (60) 

Then the solution to the equation (48) with initial condition (53) reads 

q(s) = w(s) + w(s), (61) 

explicitly, 

~Ht7 drG 1 (s,r)cos(u; n T)q n (0) 
- M Jo 

-E- r^Gx^r)™^^. (62) 

n M J UJ n m n 

It can be shown that the solution to the homogeneous equation (55) with the initial 
conditions (54) is 
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and 



u(s 



5) = i^f rfrG2(S ' r)/(T) ' {S - t] (64) 
is the solution to the inhomogeneous equation (48) with the homogeneous initial conditions 



u(t) = , u{t) = , 



where Green function G 2 (si, s 2 ) is 

G 2 {si, s 2 ) = < 



Ui(s2)u2(si)-U2(s2)ui(si) jj? ^ > S 
Ul(s 2 )«2(s2)— "1(S2)«2(S2) 2 1 



otherwise. 
Hence, we get the solution to Eq. (48) with the initial conditions (54) 



q(s) = u(s) + u(s) 
= Us) ^ 



u 1 (s)) q (t) + - lis)Pit) 



1V 7 ^ v ' u^t) M 
C n ft 

+ 5^77 / drG 2 (s,T)cos(u n T)q n (0) 

1V1 J s 



Substituting (66) into (47), one obtains 



sin(u; n r)p n (0) 



q n (t) = / ds sin[cj„(t - s)] \u 2 (s) - ^fiL^s) J g(t) 

m n uj n Jo \ Ui{t) J 

/ as sm|cj„(t — s)J — 

m n uj n Jo u 



+q n {0) cos(u n t) + 



ui(t) M 
p n (0) sm(uj n t) 



1 C C /■* /■* 

■'- ^ ; in I 



C n C„ 



f ds [ drsm[u n (t - s)]G 2 (s,r) cos(u m r) q m (0) 

Js 



~T7Z^ / " s / ^smw B t-s G 2 (s,r) p m (0) . 

M V m n u n Jo Js m m u m 



sin(a; m r) 



By using (42) we immediately arrive at 

Cl ft 



A (t) = - £ f sinK(t - s)] L(s) - ^ui(a) ) , 

^m n u n Jo \ u^t) J 

1 C 2 ft u-i(s) 
m = -T7j:— ds sinMt-s)]- 1 



(65) 



(66) 



(67) 



(68) 
(69) 
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Furthermore, A, B can be written as 

A(t) = 2 C dsrj(t - s)u 2 (s) - 2^& C dsrj(t - s)ui(s) , (70) 

JO Ui(t) Jo 

B ^=Mmfo ds,]{t - s)Uils) - (71) 

From (62), the momentum of the Brownian particle is then 

p(s) = Mq(s) 

.(, l(8) _M| i2(s) ) M9(0) + M£| p(0) 

-V / dTG[{s,T) cos(cj„t) q n (0) 

n JO 

-Y I* dTG[(s,r f nM ^ . (72) 

Here "prime" stands for derivative with respect to the first variable of Gi(s,t). 
With these results, It can be easily shown that 

C(t) = A /* d\Gi(t, X)u{t - A) 



M Jo 

2% ft 



M 2 



J* ds J* dr J* d\rj(t - s)Gi(f, \)G 2 (s,t)v(t - A) , (73) 



and 



£>(t) = ft f dXGUt, \)u(t - A) 
Jo 

-— / / dr / dWi - s)G' 1 (t, X)G 2 (s, rW(r - A) . (74) 
M Jo Js Jo 

where v(s) is defined as 

r+oo 1 

z/(s) = / dujliuj) coth.(-Tiujj3) cos(a;s) . (75) 
Jo 2 

It is seen that the coefficients A(t), B(t),C(t), D(t) are dependent only on the kernels rj(s) 
and v(s) and the initial state of the bath, not dependent on the initial state of the system. 
Once the spectral density of the environment is given, the elementary functions Ui(i = 1, 2) 
can be solved from equation (55), the Green functions Gi(i = 1,2) are then obtained. Thus 
the coefficients of master equation can be determined. 
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IV. PARTICULAR CASES 



In this section, we will consider some special cases. Let us at first treat a special case 
in which we assume that the interaction between the system and environment is weak, so 
the C n are small. In this case, the coefficients are of simple forms, and the determination 
of these coefficients is very simple and straightforward. We shall work out these coefficients 
directly using the method in the last section, rather than the general formulae. 

The solution to Eq. (46) may be written as 

p„(0) sin(w n t) 



q n (t) = q n (0) cos(cj„t) + 

m n LO r , 

C n ft .sm[w n (t - t' )] 



_^ r dtf ^n{t-t n 

m n Jo u n 

C n ft ,.am[w n (t - t')] sm[Q(t' - t)] p(t) , , „, 2 ( 



dr _j^l -^ + 0{C 2 n ). (76) 

m n Jo u n Q M v ' 

Using Eq. (37) and ignoring terms with higher than the second order of C n , we get 

Ec»<?»(f)> = {E-— f dt ' sinMt ~ r)] cosMt'-t)]}( q (t)) 

M „ m n Jo uj n is 

Then we immediately get 

Ait) = 2 f dsr](s) cos(fis) , (78) 
Jo 

m = -jkJo dsr]{s)sin{ns) - (79) 

We next evaluate J2 n C n {q(t)q n (t)) and C n (p(t)q n (t)) . After a few manipulations, we 
arrive at the expressions 

C(t) = -Y,Cn{(q(t)q n (0))cos(u; n t) + ( q ( t )p n (0)) S -^^} , (80) 
„ m n uj n 



and 



D(t) = -Y,C n {(p(t)q n (0))cos(u n t) + (p(t)p n (0))^^} . (81) 

n m n uj n 
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To calculate C(t) and D(t), we need to expand q(t) up to the second order of C n : 

,( M )cos W + ^^ 
C n ft ^ sm[Q(t - s)] 



~ 2- m Jo n cos K s ) 9n(o) 

_ \ " Q± [ 
V M Jo 



C n ft ^sin[fi(£ - s)} sin(u n s) p n (0) 



+ 0(C 2 n ) . (82) 
The expansion of p(t) is easily obtained from that of q(t) : 

p(t) = Mq(t). (83) 
With these results it is easy to compute C(t) and D(i): 

C{t) = ^uIo dSU{s)Mnsl (84) 

D(t) = nC dsu(s) cos(fis). (85) 
Jo 

This simple example exhibits the time dependency of the coefficients of the master equation 
in a general environment. Eqs. (78), (79), (84), (85) are in agreement with Hu, Paz and 
Zhang [10]. 

As another example, we briefly discuss the purly Ohmic case in the Fokker-Planck limit 
(a particular high temperature limit), which has been extensively discussed in the literature 
[6,10]. In this case one has 

r](s - s') = M-f5'(s - s') (86) 
2M 7 k B T u 

v(s — s ) = o(s — s ) (87) 

Then the equation (55) reduces to 

u{s) + n 2 Icn u(s) + 7 «(s) = -2 7 <?(s)u(0) (88) 

where Q 2 en = il 2 — 275(0). After solving this equation, a few calculations give 

15 



A(t) = -2M 7 5(0), (89) 

B(t) = 2 7 , (90) 

C(t) = 0, (91) 

D(t) = 2M 1 k B T, (92) 

then the Wigner equation reads: 

m w = ar + MfiL9 *r + 27 + 2M74bT t^ • (93) 

In this regime, the coefficients of this Wigner equation are constants. 
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